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FOREWORD 


This  report  describes  a simulation  program  and  related  verification 
hardware  developed  for  the  Air  Force  Avionics  Laboratory  by  the  IBM 
Corporation,  Huntsville,  Alabama.  This  six  month  effort  was  performed 
under  Contract  F33615-75-C-1133  and  Work  Unit  20520303.  MIL-STD-1553A 
data  networks  are  simulated  using  iterative  time  domain  calculations 
and  results  verified  with  breadboard  hardware.  Correlation  of 
simulation  results  with  hardware  is  excellent  and  modelling  Improvements 
are  presented  for  skin  effects,  transformer  coupling  and  timing. 

The  results  siininarlzed  in  this  report  were  obtained  by  Bud  Balliet, 
Program  Director,  Jame.-i  W.  Dalton,  William  W.  Scott  and  James  R.  Farl. 
Michael  J.  O’Connor  was  the  contract  monitor  for  the  Air  Force  Avionics 
Laboratory. 


TABLE  OF  CONTENTS 


Section 

Title 

Page 

I 

Introduction 

1 

II 

Theory 

4 

A.  Trenemlsslon  Line  Siauletion 

4 

B.  Reflections 

14 

C.  Stub  SiBulstion 

16 

D.  Trsnsfomer  Modeling 

19 

E.  Suaaary  of  Theoretical  Approach  and  Assunptions 

28 

III 

Prograa 

30 

A.  Organization 

30 

B.  Basic  Routines 

30 

1.  Main  Routine 

30 

2.  Subroutine  Filter 

30 

3.  Subroutine  Plotter 

30 

4.  Subroutine  SAVPNT 

41 

S.  Subroutine  PL4020 

41 

IV 

Slnulatlon  Operation 

47 

A.  DEC-10  Compability 

47 

B.  Input  Data 

47 

1.  Input  Definitions 

48 

2 . Array  Inputs 

50 

C.  Constraints  of  the  Simulation 

50 

D.  Coupler/Stub  Calibration  Discussion 

51 

V 

Laboratory  Evaluation 

53 

A.  Data  Bus  Breadboard 

53 

B.  Test  Data 

58 

VI 

Conclusions  and  Recommendations 

69 

References 

70 

iv 


LIST  OF  ILLUSTRATIONS 


Figure 

Title 

Pate 

1 

SlaulaCor  Elaacnta 

2 

2 

Tranafoner  Coupler  Schcaatic 

3 

3 

Poaalble  Equivalent  Circuit  of  a Tranaaiaalon  Line 

5 

4 

Line  Characterietlca 

8 

5 

Reaponae  of  Saaple  Cable  to  a Unit  Step  Input 

9 

6 

Line  Reaponae  to  Step  Input 

10 

7 

Skin  Effect  Line  Siaiilation  Algoritha  (I|moring 

Tl»e  Delay) 

12 

8 

Referencea  for  Eatiaating  Function 

13 

9 

Reflection  Slgnala  at  Bua/Coupler  Junction 

13 

10 

Average  Energy  Spectral  Denalty  for  MlL-STD-1553 

Syaetrlcal  Waveforaa 

18 

11 

Output  Reaponae  to  Tranaforaer  11 

21 

12 

Tranaforaer  and  Stub  Line  Equivalent  Circuit 

22Figure  13. 

13 

Equivalent  Circuit  froa  Stub  to  Main  Bua 

24 

14 

TRANS2  Teat  Setup 

26 

15 

Output  Reaponae  to  TRAMS2 

26 

16 

Block  Dlagraa  of  Data  Bua  Cable  Siaulation 

31 

17a 

Main  Subroutine  (Sheet  one  of  two) 

32 

17b 

Main  Subroutine  (Sheet  t%K>  of  two) 

33 

18 

Filter  Suhroutinea 

34 

19 

Declaal  Codea 

37 

20 

Subroutine  Ploter  (Key^l) 

38 

21 

Subroutine  Ploter  (Key-2) 

39 

22 

Subroutine  Ploter  (Key-4) 

40 

23 

SAVPNT 

42 

24 

PU020  Flow 

44 

25 

PLTEND  Flow 

45 

26 

PLOTID  Flow 

46 

27 

Coupler  Box 

54 

28 

Breadboard  Driver  Schenatic 

56 

29 

Breadboard  Driver  Input  and  Output  Waveforaa 

57 

30 

Setup  A 

59 

31 

Setup  B 

59 

32 

Slaulator/Lab  Coapariaon  for  120  Feet  of  Cable  into 

68  ohaa 

60 

33 

Siaulator/Lab  Coapariaon  for  one  10  Foot  Stub 

61 

34 

Slaulator/Lab  Coapariaon  for  one  20  Foot  Stub 

62 

35 

Slaulator/Lab  Coapariaon  for  eight  10  Foot  Stuba 

63 

36 

Slaulator/Lab  Coapariaon  for  eight  10  Foot  Stuba 

64 

37 

Siaulator/Lab  Coapariaon  Driving  Stub  #1 

65 

38 

Slaulator/Lab  Coapariaon  at  Stuba  #1  and  #8 

66 

39 

Siaulator/Lab  Coapariaon  for  one  20  Foot  Stub 

67 

40 

Siaulator/Lab  Coapariaon  for  eight  10  Foot  Stuba 

68 

V 


LIST  OF  TABLES 


Table  Title 

1 Line  Algorithm  Sxjmmary 

2 Reactance  and  Corresponding  Capacitance 

presented  by  Ten  and  Twenty  Foot  Lossless 

Stubs  Terminated  in  an  Open 


Page 

15 

19 


SUMMARY 


The  widespread  application  of  baseband  data  buses  for  avionics  system  inter- 
face has  obriated  the  importance  of  computer  simulation  for  design  and 
evaluation  of  complex  cable  netwrks.  This  report  describes  a general  purpose 
digital  computer  program  that  provides  this  capability.  Programming  routines 
are  coded  in  FORTRAN  IV  for  easy  modification  and  for  operation  on  either  an 
IBM  System  370  or  a Digital  Equipment  Corporation  DEC  10  computer.  The 
simulator  is  specifically  designed  to  emulate  one  megabit  MIL  STD  1353  type 
waveforms  and  bus  configurations. 

A typical  MIL  STD  1553  data  bus  conslrts  of  a main  bus  accessed  through  shorter 
cables,  herein  called  stubs.  These  stubs  present  a capacitive  load  to  the 
bus  coupler  which  is  located  at  the  main  bus,  stub  junction.  Coupler  trans- 
formers with  their  stub  loads  are  simulated  as  second  order  systems  with 
paraaMters  derived  from  measurement  data.  Tustins  transformation  is  used 
to  produce  different  equations  that  permit  time  domain  computation.  Itera- 
tive time  domain  calculations  are  employed.  Reflection  coefficients  are  used 
to  calculate  signal  imperfections  caused  by  junction  and  termination  dis- 
continuities. Transmission  line  filter  effects  are  accommodated  by  an 
independent  algorithm  that  takes  skin  effects  into  account.  Procedures  are 
Included  for  operating  the  program  and  for  characterizing  the  line  and  the 
transformer  couplers  from  laboratory  data.  The  program  will  permit  stubs 
of  lengths  up  to  nominally  20  feet  at  numerous  points  along  the  main  bus 
trunk.  Multiple  stubs  can  also  be  simulated  so  that  the  amount  of  stub 
loading  is  not  bounded. 

To  verify  validity  of  simulation,  a breadboard  of  a typical  bus  network  was 
constructed.  This  breadboard  permitted  a worst  case  configuration  of  a 300 
foot  main  bus  trunk  with  8 stubs  connected  via  either  of  two  types  of  trans- 
former couplers.  Results  produced  by  the  simulation  agreed  favorably  with 
data  taken  on  the  breadboard.  Additional  transfomer/stub  characterization 
study  is  recommended  to  better  correlate  transformer  parameters  to  the 
software  model. 


vil 


I. 


INTRODUCTION 


The  class  of  data  bus  defined  by  MIL-STD-1553  (reference  1)  employs  time 
division  multiplexing  at  baseband  with  a unique  synchronization  waveform 
followed  by  Biphase  L (Manchester  II)  coded  binary  data.  Tlie  main  bus 
trunk  can  be  accessed  by  as  many  as  33  cable  stubs.  A data  bus  system 
thus  represents  a complex  network  for  simulation. 

Network  simulation  can  generally  be  approached  in  either  the  time  domain 
or  frequency  domain  with  each  approach  providing  advantages  and  disadvan- 
tages. In  general,  time  domain  simulation  will  simplify  handling  complex 
waveforms  and  was  the  approach  IBM  selected  during  a Company- funded  devel- 
opment in  1974.  This  development  resulted  in  a FORTRAN  coded  simulation 
tool  that  emulated  MIL-STD-1S53  networks,  but  did  not  provide  transformer 
coupling. 

Recognizing  the  importance  of  a network  simulation  tool  for  continued  MIL- 
STD-1353  related  network  study,  this  contract  was  awarded  to  IBM  by  the  Air 
Force  Avionics  Laboratory  for  supplying  a modified  version  of  IBM's  original 
simulator.  Modifications  were  to  incorporate  transformer  coupling  and  add 
provisions  for  transmission  line  skin  effects  and  other  features  to  general- 
ly improve  utilization  flexibility. 

Figure  1 shows  the  elements  of  the  data  bus  network  simulator.  The  cable 
of  specific  Interest  is  twisted  shielded  wire  pair  and  is  shown  divided  into 
line  segments.  The  main  trunk  contains  up  to  30  segments  each  of  identical 
length.  For  a maximum  MlL-STD-1553  configuration,  each  line  segment  would 
be  10  feet  long  to  simulate  a 300  foot  main  trunk.  At  each  line  segment 
junction,  stubs  can  optionally  be  connected  through  couplers.  A coupler 
consists  of  a pair  of  isolation  resistors  and  a transformer  connected  as 
shown  in  Figure  2 . Again  referring  to  Figure  1 , each  stub  in  turn  also 

contains  cables  that  terminate  in  a load  R^.  R^  simulates  the  impedance  of 

a receiver/transmitter.  The  ends  of  the  main  trunk  are  also  terminated  in 
an  Impedance  R^,  normally  but  not  necessarily  equal  to  the  characteristic 
Impedance  of  the  line. 

Normally,  one  of  the  receiver/transmitters  on  a stub  will  be  activated  in 
the  transmit  mode.  Thus,  a generator,  G,  can  be  simulated  to  drive  the  bus 
from  one  of  the  bus  stubs.  It  can  also  be  connected  to  one  of  the  ends  of 

the  main  trunk.  The  generator  produces  a variable  three-level  format  in- 

cluding a MIL-STD-1553  waveform  to  drive  the  bus  from  an  off  state  (level 
1)  to  a plus  state  (level  2)  and/or  a minus  state  (level  3).  The  source 
Impedance  of  the  generator  is  set  at  a value  that  can  differ  from  R,  or 

Z . 

o 

As  noted  in  Figure  1 , "N"  «tubs  can  be  connected  to  a given  junction  of 
the  main  bus.  "N"  can  vary  and  can  be  any  positive  real  integer.  Two 
restrictions  are  placed  on  the  use  of  multiple  stubs  at  any  junction. 

These  are  that:  1)  "N"  must  equal  one  if  the  generator  is  active 

through  that  specific  junction,  and  2)  all  "N"  stubs  must  be  of  the 
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N = 1 . 2, 3 STUBS  AT  EACH  AiNC^’ON 
G ° GENERATOR 
Rj  > TERMINATION  RESISTANCE 
R,^  • LOAD  RESISTANCE 

Figure  1.  SlBulaCor  Eleaente 
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‘nMSTED  SHIELDED  PAIR  LINE 


STUB 

Figure  2.  Transforacr  Coupler  Schcaatlc 


< » MAIN  TRUNK 


saae  length  for  all  active  .iunctlons.  Thus,  the  nmber  of  stub  loads  that  can 
be  slsKilated  la  not  bounded. 

The  report  Is  organized  to  first  cover  basic  theory  of  network  simulation. 

This  Is  provided  In  Section  II  and  Includes  the  treatment  of  the  transmission 
line,  couplers  and  Junction  calculations.  Section  III  describes  the  basic 
program  organization  and  relates  the  theory  to  program  routines.  With  an 
understanding  of  the  program.  Section  IV  provides  operation  Instructions  and 
describes  operational  capability  and  limitations.  Section  V describes  the 
breadboard  used  In  the  simulator  development  and  documents  test  data.  The 
final  section  briefly  provides  conclusions  of  the  development  program. 
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II.  THEORY 


This  section  describes  the  basic  data  bus  elements  to  be  simulated  and  pro 
vldes  the  theoretical  basis  for  the  simulation. 


A.  TRANSMISSION  LINE  SIMULATION 

Each  of  the  transmission  line  segments  Identified  In  Figure  1 will  modify 
the  signals  passing  up  and  down  the  data  bus.  The  task  of  the  transmission 
line  simulation  Is  to  accurately  duplicate  the  transfer  characteristics  of 
a selected  line  via  a set  of  parameters  that  accurately  describe  It.  It  Is 
furthermore  highly  desirable  to  use  parameters  that  are  easily  obtained  by 
laboratory  measurement. 

1.  Applicable  Line  Parameters 

Transmission  line  theory  Is  well  documented  (see  References  2,  3,  and 
4 for  specific  examples)*.  Most  treatments  start  by  modeling  the  dis- 
tributed properties  of  a length,  i , by  equivalent  Rs,  L and  C as  shown 
In  Figure  3 . The  resistance.  Inductance  and  capacitance  Is  defined 
per  unit  length.  From  this  model  idealized  transfer  equations  and  key 
line  characteristics  are  readily  derived.  For  a line  free  of  reflection 
(Infinite  length  or  uniquely  terminated),  It  Is  shown  as  a function  of 
frequency,  that, 

^out  ^ e"^^ 


where Y Is  the  propagation  function. 

Y(u)  -^(R  + JwL)(G  + juC)  (2) 

In  terms  of  real  and  Imaginary  components,  ^ 

Y(u)  - a(w)  + J6(w)  (3) 

a Is  the  attenuation  term  and  Is  a function  of  frequency.  is  the  delay  or 
phase  term  and  also  frequency  dependent.  In  well  behaved  lines  (non 
dispersive)  B Is  linear  with  frequency  and  merely  represents  a constant 
time  delay,  T£,  In  the  time  domain. 


^Unless  specifically  noted  otherwise,  transmission  line  equations  used  in 
this  report  were  all  fully  developed  In  these  references. 
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The  characteristic  lapedance  of  a line,  the  Instantaneous  lapedance 
presented  by  the  line.  Is  defined  as 


0 y G + jwc 

It  can  be  shown  that  In  sioat  practical  lines,  the  values  of  R,  L,  G 
and  C are  such  chat  freqtiencles  above  a few  tens  or  hundreds  of  kilo- 
hertz, the  characteristic  lapedance  obtains  an  approxliaately  constant 
value.  The  value  of  Z Is  easily  aeasured  or  detenslned  from  published 
data.  ° 

If  the  paraawters  of  Figure  3 are  treated  as  lumped  values,  a time 
variant  transfer  function  Is  easily  derived.  However,  over  practical 
ranges  of  operating  frequency,  the  series  resistance  R and  Inductance 
L vary  with  frequencies  by  vlrtiie  of  the  skin  effect  phenomena. 


2.  Skin  Effect  Response 


Hlglngton  and  Nahman^  treated  the  R^  and  L terms  as  frequency  variables 
and  showed  that  for  a line  %/hoae  attenuation  characteristics  are  propor 


tlonal  to 


Che  output  of  the  cable  responding  to  a unit  step  in 


put  Is  accurately  described  by: 


f(t  + Tl)  - cerf /| 

where:  cerf  Is  the  complementary  error  function 

T ■ unit  delay  of  line 


(5a) 


o 


(5b) 


f^  “ reference  frequency 
a “ attenuation  at  f^  (Nepers) 

The  term  "B"  Is  easily  awasured  In  the  laboratory  or  determined  from 
published  data.  Although  work  performed  by  Wlglngcon  and  Nahnan  was 
specifically  fer  coaxial  cables,  the  twisted  shielded  pair  under  con- 
sideration will  behave  much  the  same  as  coax  over  a considerable  fre- 
quency range. 


^R.  L.  Hlglngton  and  N.  S.  Nahnan,  "Transient  Analysis  of  Coaxial  Cables 
Considering  Skin  Effect,"  Proceedings  of  the  IRE.  Feb  1957,  pp  16^-174. 
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Figure  4 show*  cheracterletlcs  of  the  twisted  shielded  pair  (TSP)  used  in 

the  verification  breadboard  of  a larger  cable  RG  108A/U.  The  region  at 

higher  frequency  where  the  attenuation,  Ot  is  linear  with  frequency  on  the  log/log 

plot  corresponds  to  the  condition  for  compliance  with  equatlon(5).  Thus,  the 

sample  cables  will  reasonably  respond  in  accordance  with  equation  (5)  at  high 

frequency. 

Figure  5 shows  the  time  response  of  the  #24  Teflon  TSP  sample  depicted  in 
Figure  4.  The  theoretical  curve  is  the  response  anticipated  by  equation  (5) 
using  a value  for  B derived  from  Figure  4.  The  fact  that  the  actual  response 
does  not  achieve  an  amplitude  equal  to  the  theoretical  is  attributed  to  the 
higher  attentuatlon  at  lower  frequency,  i.e. , deviation  from  the  linear 
y/lTf  region.  The  actual  response  tracts  the  theoretical  closely  along  the 
edge  of  the  input  rise  (region  of  high  frequency). 

3.  Transmission  Line  Algorithm 

In  the  simulation  program,  signals  passing  up  and  down  the  data  bus  are 
treated  as  incremental  step  functions.  For  a transmission  line  segment, 
an  input  step  of  amplitude  *'E"  will  produce  an  output  v(t-fT  ).  See 
Figure  6.  During  the  region  where  equation  (5)  applies, 


v(t4Tl) 


(6) 


This  expression  can  be  evaluated  over  the  time  Increments  between  t.,  t , t., 
t.,  etc.  vhere  each  Increment  is  of  identical  duration  and  corresponds  to 
tne  computation  time  of  the  program.  Equation  (6)  for  a single  step  becomes: 


v(t+T£)  - E (cerf + cerf^ 

+ c.rf4 

' t ^ 

ro 

+' 


■) 


(7) 


For  a given  length  of  line  and  constant  time  increments,  each  evaluation  will 
be  constant  and  equation  (7)  is  represented  by: 

v(t4T)l)  - E(Cj^  + Gj  + Gj  + ) (8) 

®1*  ^2*  *^°®*^*"^* 
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Figure  5.  Reeponee  of  Soplc  Cable  to  a Unit  Step  Input 


Figure  6.  tine  Response  to  Step  Input 
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An  arbitrary  signal  Into  a line  segment  can  be  treated  as  a sequence  of 
small  step  functions.  This  Is  depicted  in  Figure  7.  The  signal  f(t) 
la  shown  being  composed  of  superimposed  step  Inputs  ^E,  occurring  at  tj^, 

AE2  at  t2»AE^  at  t^,  etc.  From  Figure  7,  It  can  be  seen  that: 

''out  ■ “1^  '1  + =2  “(1-1)  + '3  i*(l-2)  * 

where:  1 > present  value 

1-1  ■ first  past  value 

1-2  • second  past  value,  etc. 

AE  ■ amplitude  of  incremental  Input  changes 

The  above  expression  Indicates  that  the  present  output  Is  a function  of  the 
previous  step  Inputs.  If  the  number  of  past  input  evaluations  were  made 
very  large,  the  evaluation  would  closely  approach  equation  (6).  Since  this 
is  not  practical  or  desirable,  an  estimating  function  Is  Introduced  to  allow 
the  output  to  approach  an  attenuated  value  of  the  Input  after  a fixed  number 
of  paat  value  calculations. 

The  estimating  function  calculation  Is  referenced  to  Figure  8.  It  consists 
of  the  implementation  of  the  differential  equation  for  allowing  the  output  to 
approach  the  Input  exponentially. 


dv 

dt 


k[V-v(t)] 


(10) 


In  the  diagram  of  Figure  9,  the  output  after  t 
final  value  at  V.,  passing  through  V.  at  the  t . 
Is  simply:  ^ 


Is  allowed  to  approach  a 
The  solution  to  the  above 


V 


Evaluating  the  above  at  V - (V^-Vj^) 


(11) 


T 


(12) 
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+ A E3G2^AE3G3 


4-  AE4G,  -»^AE4G2  + .. 


AT  TIMET, 


''out  -G,AEj  + GjAE,,..,,  + G3AE„.2,  + + V,.^  ,(A) 

PRESENT  1ST  2ND  ESTIMATING  PREVIOUS 

VALUE  PAST  PAST  FUNCTION  OUTPUT 

VALUE  VALUE  TOTAL 


Figure  7.  Skin  Effect  Line  Simulation  Algorithm  (Ignoring  Time  Delay) 
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Figure  8 


Rcfereneea  for  Eetlaatlng  Function 


Figure  9.  Reflection  Signals  et  Bus/Coupler  Junction 


3 


where  the  value  for  k 


The  complete  eatiaatlng  function  which  takes  the  form  of  equation  (10)  la: 


- ’1-3)“ 


where: 

'^(i-a) 

G7 

GA 


3rd  paat  Input 
3rd  past  output 
attenuation  term  for  line 
k factor  discussed  above 


(13) 


A transmission  line  Is  calibrated  at  some  convenient  length.  This  was 
arbitrarily  set  at  100  feet.  From  curves  similar  to  that  of  Figure  5( 
values  for  GA  and  G7  are  determined.  The  value  for  GA  Is  set  as  the  relative 
amplitude  after  250  ns;  the  value  of  G7  Is  set  as  the  relative  amplitude 
after  1500  ns. 

B.  REFLECTIONS 

As  a wa".'front  travels  down  the  main  bus  it  will  encounter  discontinues  at 
each  coupler  Junction  (and  possibly  at  the  end  terminations).  See  Figure  9. 
Part  of  the  incident  wave,  will  be  reflected  and  the  other  part,  4^1* 
will  continue  In  the  direction  of  the  original  wave.  A portion  of  the  wave, 
^Vg,  will  also  travel  down  the  stub.  The  reflection  coefficient,  , and 
transmission  coefficient,  , describe  the  amplitudes  of  these  waves. 


*”*o 
“ w*t 

o 

where  z - Instantaneous  Impedance  seen  at  a Junction 
p - 1 + V 


(lA) 


(15) 


The  line  algorithm  Is  sunsuirlzed  In  Table  1.  The  change  In  subscript  of 
Table  1,  equation  B as  compared  to  that  of  equation  A In  Figure  7 takes  the 
delay  of  a line  segment  Into  account.  For  example,  the  first  term  of  equation  A 
has  an  "1"  Index  whereas  the  first  term  of  equation  B is  indexed  "1*1".  The 
delay  of  one  line  segment  is  one  computation  cycle  as  Indicated  by  this  indexing. 


lA 


Table  1. 


LINE  ALGORITHM  SUMMARY 


I 


V 


out 


i 


E ■ Inputs 
V - output! 

1 ■ present  tine 


“(l-l)  * °2  “(1-2)  * ’^3  “(1-3) 


* ''>-3’  “ ’(l-l). 


(B) 


Where: 

- cerf 
Gj  - Ctf.^  -G^ 

Gj  - eetfi^f  - e«l  /j|j 

G7  • final  relative  asplitude  for  longest  pulse, 
detemined  at  value  after  1500  ns. 

Ti  chosen  to  provide  desired  amplitude 
’ after  250  ns. 


(C) 

(D) 

(E) 

(F) 
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The  Instantaneous  Inpedance  seen  at  a junction  is  the  parallel  conblnation 
of  the  line,  Zo,  and  the  stub  load. 


Thus, 


The  instantaneous  signal  anplitudes  are  as  follows: 


(16) 


AV„  - 

(17) 

R 

in 

AV^  - 

AV, 

in 

(1+^) 

(18) 

AV  - 

s 

AV, 

in 

(1+^)  f (TRANS  1) 

(19) 

trtiere  (TRAN  1}  is  the  transfer  function  of  the  coupler  as  viewed  from 
the  bus  side. 

It  is  noted  that  z is  a function  of  tine  making  many  of  the  dependent 
variables  also  functions  of  time.  This  point  is  discussed  later  In  the 
report. 

In  the  simulation  program  the  signal  wavefronts,  and  AV^  are  calculated 

and  their  summations  accumulated  each  cable  junction.  These  se 
summations  are  used  to  reconstruct  a composite  waveform  for  display. 

C.  STUB  SLMULATION 

There  were  two  approaches  considered  for  stub  simulation.  The  first,  which 
was  employed  during  the  original  simulator  development,  is  to  calculate 
wavefront  amplitudes  coning  out  of  the  transformer  coupler  and  permits  them 
tc  travel  to  the  stub  termination  through  cable  segments.  Reflection'  at 
a stub  termination  would  result  in  a portion  of  the  wavefront  travelling 
back  and  forth  on  the  stub  until  a steady  state  condition  is  reached.  Tlie 
first  approach  would  keep  track  of  these  wavefronts.  The  second  approach  is 
Co  treat  the  stub  as  an  impedance  or  load  to  the  transformer.  The  load  would 
vary  with  the  stub  length  and  is  a practical  consideration  if  load  variation 
is  not  significant  over  the  frequency  spectrum  of  interest.  This  second 
approach  has  the  advantages  of  (1) simplifying  simulation  of  any  reasonable  stub 
lengths  and  (2)  reducing  program  complexity  and  operating  cost  (run) 
time  and  main  score) . 


I 
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A typical  frequency  spectrum  of  synmetrlcal  MIL-STD-1553  waveform  (sychro- 
nizatlon  signal  plus  Bl-phase  data)  Is  shown  in  Figure  10.  A null  at 
two  megahertz  can  be  expected  with  ideal  syonetrlcal  waveforms.  In 
practice  deviation  from  the  ideal  will  result  in  even  harmonics  and  some 
energy  around  two  megahertz.  (The  fundamental  bi-phase  signalling  fre- 
quency is  one  megahertz.)  However,  it  clear  form  Figure  10  that  most  of 
the  energy  will  be  distributed  between  about  100  KHz  and  2 MHz. 

The  input  impedance  at  a given  frequency  of  a stub  as  seen  from  a transformer 
coupler  is  given  by 

*in  *0  [^z^  + tanh  yi- 


(20) 


where  ■>  the  load  impedance  presented  by  a remote  terminal. 


Since  the  remote  terminal  Impedance  is  relatively  high  and  losses  on  the 
short  stubs  of  interest  (less  than  20  feet)  are  small,  it  is  of  Interest 
to  examine  the  case  of  a lossless  line  with  approaching  infinity.  For 
a lossless  line  equation  (20)  becomes 


z 


in 


+ J z^  tan  (2r/X^)I 
+ J Zj^  tan  (2Tr/Xj^5t 


(21) 


where  X^  ■ the  wavelength  at  the  frequency  of  interest. 
Letting  Z^  approach  infinity,  equation  (21)  further  reduces  to 


'in 


-J 


tan 


^ I 

Xu) 


(22) 


Equation  22  suggests  that  for  low  loss  short  stubs  terminated  in  a high 
Impedance,  the  impedance  seen  by  the  transformer  will  be  capacitive^  inversely 
proportional  to  the  tangent  of  the  angle  determined  by  the  stub  length  and 
wavelength  of  interest.  For  small  angles  (short  stubs)  a linear  relation- 
ship is  indicated. 


Table  2 shows  values  of  impedance  and  capacitance  calculated  from  equation 
22  for  stubs  of  10  feet  and  20  feet  respectively.  Indications  are  that  an 
open  stub  can  be  reasonably  simulated  by  a capacitor  of  size  directly 
proportioned  to  the  stub  length  for  stub  lengths  less  than  20  feet.  Table 
suggests  a value  of  capacitance  equal  to  24  9,  pfd  where  Z equals  the  stub 
length. 


17 


Curve  taken 
from  Reference  6 


FREQUENCY  IIN  MHz) 


Figure  10,  Average  Energy  Spectral  Density  for 
MIL-STD-1553  Syoietrlcal  Waveforms 
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Table  2.  REACTANCE  AND  CORRESPONDING  CAPACITANCE  PRESENTED  BY  TEN  AND 
TWENTY  FOOT  LOSSLESS  STUBS  TERMINATED  IN  AN  OPEN. 

I ■ 10  feet A - 20  feat 


Frequency 

Xc  (ohms) 

C (pfd) 

50  KHz 

-j 

13,500 

200 

-j  6760 

470 

100  KHz 

-j 

6,760 

200 

-J  3380 

470 

500  KHz 

“j 

11,350 

235 

-J  674 

472 

1 MHz 

-j 

674 

236 

-J  334 

477 

3 MHz 

-j 

219 

243 

-j  99 

537 

Obviously*  the  above  analysla  la  not  rigorous  and  does  not  show  the  effects 
of  lines  with  loss  or  with  terwlnatlons  cf  finite  lapedance.  The  llatlted 
program  scope  precluded  further  Investigation  In  this  area.  Additional 
work  Is  recoiMended . 

Section  V presents  emperlcal  <Uta  that  verifies  the  accuracy  of  simulating 
a stub  as  a capacitor. 

D.  TRANSFORMER  MODELING 


A variety  of  transformer  models  are  possible  (See  References  7 and  8)* 
with  the  selection  generally  determined  by  the  nature  of  the  problem. 
Transformer  models  were  Investigated  but  an  accurate  correlation  of 
measurable  parameters  was  not  achieved  due  to  limited  contract  scope. 

The  approach  taken  was  to  model  the  transfer  function  of  the  coi^lned 
transformer  with  Its  stub  load  from  generalised  parameters.  Looking 
from  the  bus  to  the  stub,  this  was  shown  In  the  laboratory  to  be  closely 
approximated  by  a second  order  system  of  the  form. 


Gj(s) 


u g(Galn  1) 

S^  + 2Cu)  S + w ^ 
n n 


(23) 


where  s " complex  frequency 

u • iindamped  natural  frequency 
C ■ damping  coefficient 

Gain  1 - constant  determined  experimentally. 
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A typical  output  of  a tranafoner  with  a capacitive  stub  load  la  shown 
In  Figure  H>  A transfomer  looking  from  the  bus  Into  a st\ib  load  Is 

defined  as  TRANS  1.  Although  satisfactory  correlation  to  transformer 
parameters  was  not  achieved.  Figure  12  provides  a high  frequency  model 
that  produces  a transfer  function  of  the  form  to  equation  (23)  and 
suggests  that  additional  chcracterlzatlon  efforts  would  be  advantageous. 
Parameters  of  Figure  12  are  as  follows: 


R “ series  isolation  resistor 

X 

R^  ■ primary  resistance 

L >■  leakage  Inductance 

n > turns  ratio  (primary  to  secondary) 

2 

R n > secondary  resistance  as  seen  from  primary 
* 2 

C/n  - capacitive  load  of  stub  as  seen  by  primary 


The  transfer  function  Is  easily  shown 


G(s)- 


to  be. 


(24) 


where  R_*R  + R + Rn^ 

T X p s 

C ■ effective  capacitance 

The  capacitance  term,  C,  and  inductance  L,  can  be  shown  to  be  primarily 
contributed  by  the  capacitance  load  of  the  stub  and  the  leakage  Inductance 
of  the  transformer  respectively. 

Again,  referring  to  equation  23,  the  transformer  to  stub  characteristics 
can  be  adequately  modeled  from  a single  set  of  test  data  of  the  form 
shown  In  Figure  11.  The  voltage  gain  Is, 

V 

GAIN  !**’■  amplitude  at  steady  state  _ 2 

Input  amplitude  at  steady  state 

GAINl* adjusts  for  turns  ratio  and  coupler  losses. 

The  frequency  related  terms  are: 

T - periods  In  seconds 
f ■ 1/1  frequency  in  Hz 

0)  - 2IIf 


* GAIN  1 here  Is  not  identically  the  same  as  GAIN  1 in  the  program  as 

discussed  In  Section  IV. 
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Figure  12.  Transformer  and  Stub  Line  Equivalent  Circuit 
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The  dotplng  calculation  la  dcteralncd  by  knowing  the  tine  required  for  the 
aaplltude  of  the  tranalent  to  decreaae  by  50  percent  and  the  frequency 
oscillation  la  detenalned  by  the  oscillation  period.  Thus, 

e‘“  - 0.5  (25) 

where:  t “ tljse  for  transient  to  decrease  by  50  percent 


a.  (daiq>lng  constant) 

at-  0.6933 
o - 0.6933/t 

The  undamped  natural  frequency  Is: 


0) 


n 


The  damping  ratio  Is: 


(26) 


t-  J-  (27) 

n 

Since  t and  u can  be  determined  from  test  data,  the  undamped  natural 
frequency  and  damping  coefficient  can  be  calculated  froai  equations  (26) 
and  (27). 

As  discussed  previously,  signals  on  the  stub  will  be  reflected  at  the  end 
termination  and  sent  back  toward  the  main  bus  where  they  will  again  be 
partially  reflected  and  partially  transmitted  onto  the  awln  trunk.  In  the 
simulation,  the  transfer  of  energy  back  to  the  main  bua  Is  treated  as  a 
generator  driving  the  transformer  (defined  as  TRAMS  2)  In  the  opposite 
direction. 

The  transformer  model  from  the  stub  to  the  main  bus  Is  similar  to  the 
case  for  TRANS  1 except  in  this  Instance  the  load  Is  always  resistive 
and  eqtial  to  Zo/2.  The  main  bus  looks  like  two  transmission  lines  that 
are  effectively  driven  In  parallel.  A model  for  this,  the  circuit  of 
which  Is  shown  In  Figure  13,  Is  first  order.  Parameters  are  as  follows: 

- equivalent  source  Impedance 
R^  ■ Isolation  resistor  between  transformer  and  main  bus 
R^  - ZL  + R + R + R /n^  + Z/2*n^ 

c p S X 

R - transform  primary  resistance 
P 2 

L - Lp  + L /n  total  leakage  Inductance 
“ 2 * 

R^  - R^/n  transformer  secondary  resistance 
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The  solution  for  this  circuit  Is: 


G(s) 


V 

E 


s /2 
o n 

Ls  + 


2 


(28) 


This  sgsln  Is  of  the  form: 


62(8) 


(Geln  2HRoot  2) 
s + (Root  2) 


(29) 


GAIN  2 adjusts  for  the  source  Impedance  presented  by  the  stub,  the  turns 
ratio  and  miscellaneous  losses.  GAIN  2 Is  an  Important  parameter  becauae 
It  determines  the  aatount  of  energy  reflected  back  onto  the  bus  and  conse- 
quently reflects  signal  loss  due  to  the  coupling.  A method  of  determining 
GAIN  2 la  given  In  paragraph  IV-B. 

The  transfer  function  for  driving  the  data  btis  from  the  generator  located 
at  a stub  end  Is  Identical  to  the  model  described  above  except  for  parameter 
changes.  Parameter  changes  result  because  of  a different  signal  source 
Impedance.  This  transfer  function  Is  again  first  order  and  given  as 
follows : 


G3(s) 


(Gain  3) (Pole) 
s + (Pole) 


(30 


where  POLE  Is  related  to  Rt/L 

Using  the  test  setup  shown  In  Figure  14,  the  gain  (GAIN  3)  and  POLE 
(R  /L)  are  easily  determined.  Referring  to  Figure  15,  these  parameters 
arc  as  follows: 

1 . output  amplitude 
Input  amplitude 

T^^  - Rise  time  (63%  of  final  value) 
p.i.  - ’ 't/'' 


A bilinear  Z - transformation  Is  used  to  convert  transfer  functions  G(s) 
to  difference  equations  used  for  time  domain  calculations.  Tustln  trans- 
formation (Reference  8)  was  used  successfully  on  the  previous  simulator 
development  and  was  employed  for  transformer  transfer  calculations.  The 
Tustln  transformation  requires  substitution  for  the  complex  frequency 
variable  s. 


8 


8-1 
8 + 1 


] 
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\ 

4 


I 


vhere  T !■  the  computation  cycle  time 

With  this  substitution,  equation  (23)  becomes. 


C1F3*2  + ClF2*s  ClFl 


where:  ClFl' 

C2F2' 

C2F3 

C2F2' 

V 

E 


w2*t2-2  w*T+l 
2 2 

2*W'*T  -2. 

1.+2 

GAINl*WeT^ 

output 

Input 


Performing  cross  multiplication  yields: 

(C1F3*Z^+C1F1*  *Z4C1F1  ’ ) V-C1F4  ’ (Z^+2Z+1)  E 

Dividing  by  C1F2*Z^  and  solving  for  V yields 

V - C1F4*E+2*C1F4*Z”^*E+C1F2*Z“^*E 
-C1F2*Z"^*V-C1F1*Z"^*V 

Transforming  the  above  equation  to  the  discrete  time  domain  where  Z°  cor 
responds  to  the  present  value,  Z~1  corresponds  to  first  past  value  and  Z 
corresponds  to  second  past  value  yields 

- ClF2*[E^+2E^_j^+E^_2l"‘^l^*V^'^ 

-C1F1*V^_2 


Making  the  Tustln  transformation  Into  equation  29  yields: 


. V . C2F3(z-H)_ 

E C2F2*z  + C2F1 


where:  C2F1  - R00T2*T-1. 

C2F2  - l.-HUX)T2*T 
C2F3  - GAIN2*ROOT2*T 
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Cross  nultlpllcstlon  snd  dividing  by  C2F2*Z  yields: 

V - C2F3*  (1+*"^)E-  C2Fl*s“^ 

Now,  %nrltlng  the  sbove  equation  In  the  discrete  tine  donaln  yields: 

- C2F3*  [E^  + E^_^]  - C2Fl*Vj_j^ 

E.  SUMMARY  OF  THEORETICAL  APPROACH  AND  ASSUMPTIONS 

Trsnsnlsslon  line  elements  of  a data  bus  network  are  simulated  by  a 
special  algorithm  that  takes  skin  effect  Into  account.  All  segments  of 
the  main  bus  are  of  Identical  length,  ten  feet  for  MIL-STD-ISSS  network. 

The  algorithm,  summarized  In  Table  1,  Introduces  no  Inherent  delay. 

The  delay  of  a transmission  line  segment  becimes  the  assigned  value  for  a 
computation  cycle  and  Is  set  equal  to  the  group  delay  characteristics  of  the 
line. 

An  Input  to  a line  segment  produces  an  output,  which  during  the  next 
computation  cycle,  becomes  the  Input  to  the  next  line  eleisent. 

When  a transformer  coupler  Is  attached  to  a Junction  on  the  suiln  bus,  a 
discontinuity  Is  Introduced.  This  discontinuity  Is  slsnilated  by  a reflection 
coefficient  as  defined  by  equation  14. 

Stubs  are  normally  terminated  In  high  impedances.  Stubs  of  relatively  short 
length,  20  feet  or  less,  will  reflect  a capacitive  load. 

At  each  line  element  along  the  laaln  bus,  a transfonser  coupler  can  optionally 
be  attached.  All  transformer  couplers  are  Identical.  The  transfer  function 
of  a transformer  coupler  working  Into  a capacitive  stub  load  Is  simulated  as 
a second  order  system,  the  paraiseters  of  which  are  determined  by  standard 
techniques  from  laboratory  measurement.  This  approach  adjusts  for  transformer 
parameters  that  are  difficult  to  measure  and  the  fact  that  stubs  have  finite 
termination  and  loss  characteristics. 

Basic  simplifying  assumptions  for  the  simulation  are  as  follows: 

1.  The  transmission  line  Is  non-dlsperslve.  This  assusq>tlon 
appears  reasonable  because  dispersive  lines  would  mate 

poor  data  bus  candidates.  Most  transmission  lines  reasonably 
meet  this  assumption. 

2.  The  characteristic  Impedance  of  the  line  Is  constant  and 
resistive  over  the  spectrum  of  Interest.  The  line  used  in  the 
breadboard  did  not  fully  meet  this  criteria,  but  adverse 
effects  were  minimal. 
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3.  Stubs  present  s predoalnsntly  cspscltlve  load.  This  sssuaptlon 
should  be  valid  for  esses  under  twenty  feet  where  stub  terwinstion 
inpedsnee  is  high.  MIL-STD-1SS3  tersinstes  dictate  a high  inpe- 
dance . 

4.  The  second  order  systea  sisiulation  of  a coupler  with  its  stub  load 
accurately  siaulates  the  transforaers  under  investigation.  This 

later  assuaption  has  proven  reasonably  accurate  for  the  two  transforaers 
tested,  but  should  be  verified  for  other  saaples.  The  FORTRAN  IV  coding 
used  in  the  program  makes  alternate  transformer  models  easily  inserted. 
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III.  PROGRAM 


A.  ORGANIZATION 

The  progrea  la  structured  to  Independently  ■alntaln  the  tiae/aaplltude 
phasing,  during  aultlpsth  signal  propagation,  of  data  transalsalons  and 
resulting  reflections.  Figure  16  shows  the  signal  functional  flow. 

Incident  data  wavefronts  propagate  from  left  to  right  In  filter  sequence 
FILl  through  FIL30.  FIL  functions  lapleaent  the  transalsslon  line  segaents. 
Transalsslon  reflections  travel  froa  rl^t  to  left  in  descending  order 
through  FIL60  to  FIL31.  Twenty-nine  ISTUB  filter  paths  corresponding  to 
the  data  bus  stubs  connect  the  transalsslon  sequence  and  the  reflection 
sequence.  Data  traveling  down  the  ISTUB  paths  change  froa  the  transalsslon 
aode  to  the  reflection  aode  atreaa  at  the  aid  point  (Labeled  End  of  Bus) . 

A typical  stub  path  starts  through  TRANSl,  the  transfer  function  of  a coupler 
froa  bus  to  stiA.  TRANSl  provides  the  capacitive  loading  of  the  stub.  The 
output  of  TRANSl  travels  through  a FIL  section  to  the  stub  teralnatlon 
Indicated  by  the  FA  calculation.  F4  la  the  stub  end  reflection  coefficient 
which  deteralnes  the  portion  of  the  signal  to  be  sent  back  toward  the  aaln 
bus.  This  reflection  la  passed  through  a FIL  eleaent  and  through  TRANS2 
which  provides  transfer  characteristics  froa  stub  to  bus. 

The  prograa  Is  organized  to  aalntaln  signal  "book-keeping"  for  each  data 
bus  location  versus  tlae  and  to  call  block  function  subroutines  In  proper 
sequence.  Data  transfer  between  routines  Is  accoapllshed  through  C(MM0N. 

B.  BASIC  ROUTINES 

Three  types  of  routines  are  basic  in  the  prograa  organization.  The  MAIN 
routine  perforas  calculations,  controls  the  sequencing  of  the  calculation 
subroutine  FILTER,  and  also  calls  output  subroutines  PLOTER  and  SAVFNT. 

1.  Main  Routine 

The  book-keeping  function  associated  with  aalntalnlng  Initial  and 
subsequent  values  of  prograa  variables  and  of  sequencing  subroutines 
Is  accoapllshed  In  "MAIN."  In-line  calculations  arc  also  perfonMd 
to  generate  the  transaltted  data  string.  See  Figure  17  (2  sheets) 

2.  Subroutine  Filter 

Subroutine  Filter  perforaed  the  difference  equation  calculations 
necessary  to  advance  tlae  (T).  See  Figure  18. 

3.  Subroutine  Ploter 


The  PLOTER  subroutine  Is  designed  to  Interface  with  the  host  systea's 
plot  package.  A PL4020  subroutine  Is  supplied  to  obtain  printer 
plots  only.  The  PL4020  routine  can  be  replaced  with  the  host  systea's 
plot  aodules.  Care  should  be  taken  that  cosaunlcatlon  variables 
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TRANSMITTER 
INPUT  OPTION 


^Typical  Stub  Generator  Entry 


r 


MAIN 


•HFAOI  INtSFCT  J 

WRITE!  I OUT  • 280  ITtCr  MAX.  7.x. LENGTH!  I ) 
K500.K1 500.  alpha. R1 30 .4) 

1 ,R! 31 . 4 ) . 7C.ZL! 1 ) ^ 

FORMAT! 'OSFCT  DA TA • . 2X . I IF  I 0 . 5 ) 
READ! I N.TRANl  ) 

WRITE! lOUT .TRANl ) 

READ! 1N.TRAN2) 

WRITE! inUT .TRAN^I 
READ! IN  .GPNXFR  I 
WRITE! ICUT.GENXFR) 

READ! I N .STUB ) 

WRITE! lOUT.STUBI 
READ! IK.GFNEH I 
WRITE! lOUT .GENFR | 


I 

S S S S CALCULATE  REFLECT  FACTORS  » A I t 
• DO  270  1 = 1 .29 

IF! ISTUO! I ) .LT.  1 )GO  TO  270 
R! 1 . I > = RFFl 

R! I .2 ) = REFl 

XK!  I ) = I STUB! I I 

RP!l)=70/2. 

R!I.4)=  ! ZLin  - Z0>/!2L!I)  ♦ ZO  » 

270  CONTINUE 


**««*«*•  TRANSFORMER  DATA 

T = T/2. 

DO  5b  I I = l . 2 9 

IFII5TUBlII.LT.il  GO  TO  559 
XL  ENG  = LE  NOTH  ! 1 ) 

ODl  = WNS01*XLENG/LSTU0 ! I I 
GO  TO  560 
559  UDl  = WNSOl 
5b0  lJD2  = Twn/Wl 
003  = 1.0 

004  = DO  I*  GAIN  I 

CIF31I)  = 00  3 T*OD2  ♦ T*T*ODI 
C1F21I)  = !2.*Dr)l*T*T  - 2 . *003  I 7 C ! f J ( I ) 

CIFUI)  = !UDIAT*T  - n02*T  ♦ DD3»/CIF}(I) 
C1F41I)  = 004*T*T/CIF3! I ) 

561  CONTINUE 

*«««*«« « TRANSFORMER  DATA  FROM  STUB  TO  BUS 
C2F2  = 1 .♦TAR00T2 

C2F1  = ! T*R00T2-l . J/C2K2 

C2E3  = GAI N2*MOOT2*T/C2F2 
**  4*41  FIRST  ORDER  SIM  FOR  GEN  ErOM  STUB* 

C3E2  = I ♦ TAPOLE 
C3FI  = ITAPOLE  - II/C3F2 
C3F3  = GAIN3AT*  POLE  / C3F2 
T = T *?. 
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Figure  17a.  Jlain  Subroutine 


Sheet  1 of  2 


• r»«pi»T“v*Lce5' 


•>9  C''mtN»>F 

no  Trii  N- 1 .12^ 

f>j  lor  iM.4 

'<=1-1 

JOO  FlLC^il.ui  =F|LIM«la<) 


ICO  COMTINUE 


DLL <01 
F1L(01 
o<  L( 

F l4-<  OC 


0 1 . n 1 1 « r 1 L ( 3 1 . 2 . n -• 

?.-:{;l'.-:(F!L(7J:2.n  - 

oC  . I •!  ) = / »n&0t  ^ • I ) t_  I * * J — — — 


C ••  SISNAt  VifNfCAj 


i i ^ J 


. 1 • 

. ) . M 


t « < M/  N,  I . I I ♦ 


Figure  17b.  lUln  SubroutlM 


• •••••••••HUS  C*BLF  S!»*ULATIC*<  •••••••••••••• 

_ — - 

orLFKl)  * r ILI  inTlT^^TCri.r;?! 

O^LFPCIJ  * riL«ltlt?J  - F|LCI.l«l» 
oEL^Jin  = MLUti.n  - 

E IL(  I .2.1  I * &<  I . n •OELFH  I > * Cl  I ,2)«t5ELF2l  M 

1 ♦ 01 1 . j)»ocCr  j(  I •♦.ii  I .•  1*1  cn  .7  ••f^li  I • I .*1  - eilii.?.«h 

2 • F IL ( 1.2.21 


• •••••••••••a  THANbFCHMcH  :>|U3  SIM  rurm  riUS  TLi  STuH  •••• 

00  *25  1*1. 30_  . - 

•»'•  r*6e 

L * I ♦10 

TnSll.l.>|:Cira(|t«(FlL(M.l.t>«2.*F|L|M.I.?l»F|L(M.|,3|| 
i -cir2  1 I ) «THS1  I . 1 ,2  l-cir  11  I l«TO.;(  I . 1 . •«)  ' -■  - 

DELFKU  * T^wsd.l.ll  - TPS(t.l.2t 

neuFaiii  = THsii.i.Pi  - TRS(t.i.'>i 

0'='LF.T<n  S TOSlI.l."'!  - ThSII.I.*! 

F|LI'«.2.lt  = Cl  I . I I «iir  i FI  (I  I ♦ Cl  I .2»aD?l.F2(  1 J 

1 • Cl  I . <1  '1  1 I ♦ Cl  I . «l  aid  I .7|aTOSl  1.1.41  - FIL1M.2.4)I 

2 ♦ F ILl M.2.?  ) 

TRSI  1.1.4.)  ^ TOSH  .1.31  ' 

iHsii.i.j)  =~TPsii.r.2)  ■ 

♦ OC  T:JS1  1.1.2)  S TRSI  1 .1.1) 

42s  '"ONT  IMI. 


IF(  ISTUUIKl.EO.O.O) 


• •••  ••••  •!;.  •••SF‘*PVSP  4l»40  SfUa  S)V  FPr*M  STUt)  IT  f)US  •••••• 

O'J  62?  I»>r.TO 
K*  I’-IO 

M * I .^0 

■ ■ L 1 » 60 

FILIL.1.1)  * FU.lL.l.n  ♦ |F|LlM,2.t>  - F |L  1 M.  2 . 2 ) ) aP  I K .4  » 
OELFUI)  s FILlL.l.ll  - F|L(U.).2> 

n'-LF?n)  » FILI1..1.2I  - P1L(L.1,-»I 

OiLF.Hn  *~F1lI  lVI  .3)  -'FIHL.1.41 
FILLIL.2.11  • Cl  I . I l40Et.Fl  1 I ITGI  t .2>*/>Et.F21  I) 

1 ♦ G(  I . j)  •r>F(.F31  1 I ♦ Gl  I .4)  *1G1  1.7I4F  Il.lt.  1 .4)-F  ILLIL.2.4)  » 

2 ♦ riLtU-.2.2l  , 

HLL(t.2.f)  • FlttlL.2.3)  ■ 

F|ttlt.2.3)  4FILL1L.2.2) 

FlttlL.2.21  sF|ttlL.2.1) 


■ ••••••••a  RCFLFCT  ION  F^TOP  FCQH  STUB  tQ  6U5 

RIK.H  • PEF3 
FltPlK.I.II  zF|tLlt.2.l) 

TPSI I . 1 . I l«C2F  ia|F ILPIK. I.l )«F ttPIF. 1 .21 ) 

I -C2F1 ariSI  I . 1 .2)  

FIL“1»  . 1 . J)  =F  H.OK.1  .2J 

FILPIF.I.2I  > FILP'K.I.II 
FItIt .2.1  ) * TPSI I ...I  I 

FltlL.2.21  * T«?(I.I,2I 

FILIL.2.31  * TRSII.1.3) 

TOSlI.l. 31  4 TOSH. 1.21 
T)I.SUit.2)  ■ TRSII  .!.!> 


ri(ur«  18.  Filter  SubroutlMs. 


are  consistent  if  a change  la  made. 

a.  Internal  Parameters 

The  routine  must  be  supplied  information  in  tabular  form  on 
what  variablea  are  to  be  plotted  and/or  printed.  If  plots  are 
requested,  then  information  pertaining  to  the  plots  desired 
must  be  supplied  by  input  cards. 

No  output  parameters  are  returned  from  this  subroutine.  The 
output  of  the  routine  will  be  the  requested  plots. 

b . Input  Cards 

The  supplied  input  cards  contain  information  required  to 
obtain  a plot  of  one  or  more  program  variables,  using  the 
PLOTER  subroutine  for  the  supplied  plot  package.  In  order 
to  use  the  PLOTER  subroutine,  the  user  must  supply  a sub- 
routine named  SAVPNT.  Thia  subroutine  has  the  calling 
sequence  CALL  SAVPNT  (AREA,  2)  %rhere  the  array  AREA  can 
obtain  a maximum  of  2S  variables. 

Control  Card  (one  card) 

Columns  1-2:  Number  of  variablea  (must  ba  25  or  lass). 

For  clarity,  it  may  coincide  with  the  dimension  of  array 
AREA. 

Coluimis  3-4:  Number  of  plots  requested. 

Columns  5-6:  The  Increment  to  determine  the  number  of 

points  to  be  plotted;  e.g.,  2 will  cause  every  second  data 
point  to  be  plotted.  This  will  default  to  1 if  the  field 
is  blank. 

Columns  7-8:  The  FORTRAN  unit  number  of  the  peripheral 

storage  device.  This  will  default  to  1 if  this  field  is 
blank;  therefore,  a GO  FTOlFOOl  DD  statement  must  be 
specified  for  the  unit. 

Plot  Requests  (one  to  four  cards) 

Columns  1-2:  An  integer  value  which  denotes  a subscript 

in  the  AREA  array  containing  the  variable;  e.g.,  02  would 
denote  AREA  (2).  This  is  the  Independent  variable  for 
the  curve  (abscissa). 

Columns  3-4:  An  Integer  value  denoting  s subscript  in  the 

AREA  array  of  the  dependent  variable  to  be  used  for  the 
curve  (ordinate). 
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Column  5:  Blank 


Column  6:  The  mode  of  the  grid.  (Hode  - 0;  linear  grid 

for  both  the  Independent  and  dependent  variable.  Mode  ■ 1; 
log  grid  for  the  independent  variable  (abscissa)  and  linear 
grid  for  the  dependent  variable  (ordinate)  . Mode  ■ 2;  linear 
grid  for  the  Independent  variables  and  log  grid  for  the 
dependent  variables.  Mode  *•  3;  log  grid  for  both  the 
Independent  and  dependent  variables.) 

Column  8:  Frame  advance.  (0  - advance  to  a new  frame;  and 

1 - plot  on  the  same  grid.  To  reflect  meaningful  information 
when  using  1,  the  independent  variables  for  the  curves  should 
be  in  the  same  units.) 

Column  9-10:  A decimal  value  for  the  plot  symbol.  This  is  the 

symbol  used  to  define  the  curve.  (See  Figure  19  for  the 
decimal  values  of  the  plot  symbols.)  Unprintable  characters 
will  yield  blanks  when  using  the  printer  plots. 

The  preceding  10  character  sequence  should  be  continued  across  the  card(s) 
until  M fields  have  been  completed,  where  N ■ the  number  of  curves  requested, 
as  denoted  in  columns  3-4  of  the  Control  Card. 

Identification  Frame  Data  (one  card) 

Columns  1-72:  Contain  the  Information  which  is  to  appear  on 

the  identification  frame  of  the  plotter  output. 

Frame  Title  Data  (one  card) 

Columns  1-72:  Contain  the  heading  information  which  will 

appear  on  all  the  frames  of  the  plotter  output. 

Axis  Label  Data  (one  card  for  each  variable  - 25  maximum) 

Columns  1-72:  Contain  the  labels  of  the  coordlxute  axes. 

The  number  of  cards  in  the  set  must  be  equal  to  the  number 
of  variables  denoted  on  the  control  card.  (If  the  user 
requests  multiple  curves  per  frame,  the  axis  will  be  labeled 
with  the  label  of  the  last  variable  requested  for  that  frame.) 

c.  Flow  Charts 

The  flow  charts  of  logic  paths  of  the  subroutine  are  supplied 
on  the  following  pages.  There  are  three  logic  paths  depending 
upon  the  value  of  the  parameter  KEY.  (See  Figures  20,  21  and 
22.) 
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Figure  19.  Decimal  Codes 
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SUBROUTINE  PLOTER  (KEY-11 


Figure  20.  Sidiroutlne  Ploter  (Key-1) 
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SUtnOUTINEfLOTER  (KEY-2) 


Figure  21.  Subroutine  Ploter  (Key*2) 
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MMUTWf  riOTin  iuv-«i 


Figure  22.  Subroutine  Ploter  (Key*A) 
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4.  Subroutine  SAVPNT 

The  SAVPNT  routine  allows  the  user  flexibility  In  the  data  values 
he  wishes  to  save  for  later  plotting.  The  array  PLT  Is  filled 
with  the  variables  to  be  plotted  and  transferred  to  subroutine 
PLOTER.  Flow  charts  for  SAVPNT  are  shown  In  Figure  23. 

5 . Subroutine  PL4020 


This  routine  Is  a FORTRAN  routine  designed  to  generate  X-Y 
rectangular  graphs  on  a 125  column,  high  speed  printer. 

This  routine  can  be  replaced  with  a host  system  plot  package. 

Care  should  be  taken  to  guarantee  correct  communication  variables. 

a.  Using  Printer  Plot  Routines 

1.  Initialize 

CALL  PLOTID  (IDF,  1,  1) 

where:  IDF  • Starting  location  of  72  character  label 

for  ID  frame. 

2.  Plotting 

CALL  PU020  (NPLOT,  MODE,  NCHAR,  NP,  X,  Y,  XMIN,  XMAX, 
YMIN,  YMAX,  BCDX,  BCDY,  IDF,  lERR) 


where : 

X,  Y 

NP 

BCDX,  BCDY 
NPLOT 

MODE 


are  the  names  of  the  arrays  containing 
the  X and  Y coordinates  respectively. 

the  number  of  points  to  be  plotted. 

BCD  labels  for  the  X and  Y axes 
(maximum  of  72  characters  stored 
In  FORTRAN  "A"  array) . 

determines  the  nusiber  of  curves  per 
grid.  (NPLOT  - 1 will  advance  the 
frame;  NPLOT  2 will  use  the  previous 
frame) . 

determines  the  type  of  grid  that  will  be 
drawn.  (MODE  * 1 linear  grid  for  both  X 
and  Y;  MODE  > 2 log  grid  for  X,  linear 
grid  for  Y;  MODE  ■ 3 linear  grid  for  X, 
log  grid  for  Y;  MODE  - 4 log  grid  for 
both  X and  Y. 
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SAVPNT 


Figure  23.  SAVPNT 


42 


NCHAR 

this  Integer  selects  the  plot  symbol 
to  be  used.  (See  table  of  decimal 
codes . ) 

XMIN, 

XMAX 

mlnlanim  and  maximum  values  for  the  X 
coordinates. 

YMIN, 

YMAX 

minimum  and  maximum  values  for  the  Y 
coordinates. 

IDF 

location  of  72  characters  to  be  printed 
as  a heading  for  each  frame. 

lERR 

Thie  paragrsph  is  set  by  the  subroutine 
to  Indicate  an  error  condition: 

lERR  “ 1 NormAl  return. 
lERR  * 2 Subroutines  unable  to 

construct  readeble  grid. 
lERR  • 3 Off-scsle  plot  points  were 
encountered 

NOTE:  lERR  should  be  set  to  zero 

Initially. 

3 . Temlnating 

CAU  PLTEND  (1) 

The  plots  are  output  each  tine  a subsequent  call  is  made. 
Thus,  after  the  call  to  PLOTID  Is  msde,  there  Is  no 
output  until  the  first  call  is  made  to  PL4020  with 
NPLOT  ■ 1.  On  the  next  call  to  PLA020,  the  graph  from 
the  first  call  is  output.  On  the  call  to  PLTEND,  the 
final  graph  is  output.  The  system  of  always  seeming  to 
be  one  behind  was  necessary  to  be  able  to  have  multiple 
plots  on  a single  graph. 

b.  Flow  Charts 

Flow  diagrams  for  PLA020  subroutines  entry  points  are  shown 
In  Figure  2A,  25  and  26. 
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Figure  24.  PL4020  Flov 
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Figure  25.  PLTEND  Flow 
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IV.  SIMULATION  OPERATION 


This  section  of  the  report  provides  Instructions  for  operstlng  the  network 
slnulstor.  The  program  has  been  developed  on  an  IBM  370  for  anticipated 
use  on  AFAL's  DEC-10. 

A.  DEC-10  COMPATIBILITY 

Operational  requirements  of  the  DEC-10  are  that  they  contain  a basic  FORTRAN  IV 
compiler.  Include  peripheral  storage  of  about  5 million  bytes  and  have  a printer 
capable  of  125  characters  per  line  or  greater.  These  requirements  are  met  by 
the  model  presently  at  AFAL. 

Preliminary  tests  were  made  on  the  DEC-IC  with  program  cards  developed  for 
IBM's  original  network  simulator.  These  tests  resulted  In  minor  changes 
whlc'i  compensate  for  the  fact  that  the  DEC-10 's  FORTRAN  compiler  Is  soisewhat 
more  restrictive  than  that  used  on  the  original  development.  Rearrangesient 
of  source  cards  easily  rectified  the  minor  dlff Icultltles.  Since  the  data  bus 
network  simulator  Is  programmed  In  basic  FORTRAN,  no  program  logic  changes 
were  required  to  run  on  the  DEC-10.  Both  machine  use  EBCDIC  punch  cards  as 
their  normal  read  mode.  Disk  space  on  the  DEC-10  was  sufficient  for  the 
simulator. 

The  acceptance  test  demonstration  was  accomplished  on  September  4,  1975.  Again, 
only  minor  card  changes  were  required  to  adjust  for  operating  differences  In 
the  two  computers. 

B.  INPUT  DATA 

Program  Input  Is  accomplished  by  using  the  FORTRAN  NAMELIST  format.  The  Input 
data  shown  below  Is  sufficient  to  characterize  the  data  bus  components  and 
configuration. 

&SECT 

T-16 . E-9 , CCMAX-300,ZX-112  ,LENGTH-120*10,K500-0 . 942809467  ,K1500-0 . 954521364 
ALPHA-. 009, R(30,4)-0.0,R(31,4)-0.0,Z0-68.,ZL-31*2200. 

&END 

&TRAN1 

GAINl-0 . 25 , WNSQl-5 . 03E14  ,TW0ZW1-8E6 , IDEALT-0 

&END 

&TRAN2 

GAIN2-0.1,ROOT2-2E7 

&END 

&GENXFR 

GAIN3-0 . 38 , POLE-8 . 9E6 , REFl— . 044 , REF3-0 . 0 
&END 
&STUB 
NPSKIP-2, 

ISTUB-0, 0,0, 1,27*0, 

LSTUB-0,0,0,10,27*0 

&END 

&GENER 

GENMAX-10.  ,TSLOPE-50. E-9 ,NCCIN-94 ,NCCF-32 , IGEN-0,ZG-68 

&END 

&SKIN 

Gll- . 973 , G12- . 0079 , G13- . 0060850 
&END 
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1.  Input  Definitions 


The  program  variable  definitions  and  functional  groupings  are  defined 
below.  Unless  otherwise  Indicated,  Impedances  are  in  ohms,  and  time 
is  in  seconds. 

SECT  (This  group  defines  transmission  line  dimensions  and  parameters) 

T Sample  data  period  or  computation  cycle  time.  It  Is 

set  equal  to  the  transport  delay  of  a section  of  line 
on  the  main  bus. 

CCMAX  The  number  of  computation  cycles  to  terminate  a run. 

This  Is  generally  determined  by  the  length  of  time 
desired  for  display  consistent  with  plotter  capability. 

ZX  The  stub  isolation  resistance. 

LENGT  The  length  of  segments  modeled.  This  entry  signifies 

the  total  number  of  segments  (e.g.,  120;  see  Figure  16) 
times  the  length  of  each  in  feet  (e.g.,  10  feet). 

KSOO  The  relative  output  voltage  to  input  voltage  of  100  feet 

of  cable  terminated  in  ZO  to  a step  input  after  250  ns. 

See  point  indicated  in  Figure  5. 

FT500  The  relative  output  voltage  to  input  voltage  of  100  feet 

of  cable  after  1500  ns.  See  point  indicated  in  Figure  5. 

ALPHA  Line  attenuation  in  db  per  foot  at  one  megahertz. 

R(30,4)  The  reflection  coefficient  at  the  left  end  of  the  main  bus. 

R(31,4)  The  reflection  coefficient  at  the  right  end  of  the  main 

bus. 


ZO 

ZL 


The  characteristic  impedance  of  the  cable. 

The  impedance  of  a stub  termination.  This  is  an  array 
that  permits  different  Inputs  at  each  termination. 

(See  paragraph  IV-B-2.) 


TRANS  1 (This  group  defines  characteristics  of  couplers  as  viewed  from 
the  bus) . 

GAIN  1 A gain  factor  as  described  in  paragraph  II-D.  Because 
of  the  programming  technique  used,  the  entry  should  be 
multiplied  by  0.5/(l+Ref  1). 
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VfNSQl  The  natural  frequency  of  oscillation  squared  (see  Paragraph 
II-D.) 

TWOZWl  The  term  2 Qui  in  the  denominator  of  equation  23  (See  Paragraph 
II-U.) 

TRAN  2 (This  group  defines  characteristics  for  transferring  reflected 
energy  back  to  the  main  bus  when  a stub  does  not  contain  the  generator) 

GAIN2  The  gain  factor  as  described  in  Paragraph  II-D, 
equation  29. 

ROOT2  The  Rt/L  tern  appearing  In  equation  29. 

GENXFR  (This  group  generally  defines  transformer  characteristics  for 
transmitting  from  a stub) 

GAIN3  The  gain  factor  similar  to  GAIN2  above. 

POLE  An  R/L  term  similar  to  ROOT2  above 

REF  1 r^^iThe  reflection  coefficient  Introduced  at  each  stub 
location.  Equation  14  is  used  to  determine  this  entry. 

STUB  (Inputs  here  are  for  data  bus  and  plot  configuration  control) 

NPSKIP  This  entry  permits  fewer  plotting  points  than  normal. 

Resolution  is  com,>romlsed  for  viewing  a longer  time  base. 

The  sample  of  NPSKIP* 2 means  every  other  point  is  printed. 

ISTIIB  An  array  that  defines  the  number  of  stubs  at  each  of  the 
29  main  bus  locations.  See  Paragraph  IV-B-2  that  follows. 

An  entry  of  2 means  two  parallel  stubs  at  a given  location. 

LSTUB  An  array  that  defines  the  length  of  stubs  at  each  of  the 
29  stub  locations.  For  example,  0 Indicates  no  stubs 
and  3 indicates  3 stubs. 

GENER  (Generator  characteristics  are  in  this  group) 

GENMAX  This  entry  defines  the  output  voltage  of  the  generator. 

TSLOPE  The  rise  and  fall  times  of  the  generator  signal  as  defined 
by  a signal  from  zero  volts  to  GEIMAX. 

NCCIN  The  number  of  computation  cycles  in  each  half  of  a 
synchronization  waveform  which  always  appears  as  a 
positive  followed  by  a negative  level  of  equal  pulse 
width. 
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NCCF  The  maiber  of  computation  cycles  In  each  bl>phase  data 

pulae  that  follows  the  second  half  oa  a synchronization 
signal. 

IGEN  The  location  of  the  generator. 

0 - nain  bus,  left  end; 

1 “ stub  location  one,  2 ■ sub-location  two 
3 - stub  location  three,  etc. 

ZG  The  source  Impedance  of  the  generator. 

SKIN  (This  group  provides  entries  for  calculation  related  to  complimentary 
error  functions) 

Gll,  G12,  G13  Values  of  Gl,  G2,  and  G3  are  calculated  from  equations 
given  In  Table  1.  The  cable  parameter  "B"  Is  as 
discussed  on  page  6. 


2.  Array  Inputs 

Several  of  the  Inputs  described  above  are  arrays.  An  array  consists  of  "N" 

consecutive  Inputs  of  the  sane  type.  For  example,  ISTUB  defines  the  location 

and  number  of  stubs  at  each  of  the  29  possible  locations  along  the  main  bus. 

This  Input  can  be  Introduced  as  twenty-nine  consecutive  integers  such  as  follows 

0,0,1,0,0,1,1,0. . . (29  entries) 

An  array  can  also  be  Inttroduced  by  an  abbreviated  method  such  as: 
2,0,2,20*0,0,2,2,3*0 

Where  20*0  and  3*0  Indicate  20  and  three  consecutive  locatlona  having  no  stubs. 

C.  CONSTRAINTS  OF  THE  SIMULATKXN 

As  presently  written,  the  simulation  program  has  the  following  characteristics: 

1.  The  main  bus  Is  made  up  of  30  cable  segments  of  identical  length.  The 
length  of  a csble  segment  can  be  changed  however,  by  modifying  T. 

2.  All  couplers  must  be  identical. 

3.  The  generator  routine  produces  a NIL-STD-1553  format  consisting  of  a 
S3mc  waveform  followed  by  consecutive  bi-phase  L "ones".  Other  patterns 
%rould  require  a program  change.  If  the  signalling  spectrum  Is  modified 
considerably,  stub  simulation  techniques  should  be  reinvestigated. 

4.  The  stub  termination  Impedance  should  be  about  IS  Zo  or  greater  for 
accurate  stub  simulation. 
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5.  Th«  length  of  a stub  for  display  tine  phasing  la  always  equal  to  that 
of  a line  elenant  on  the  sain  bus.  This  Introduces  a snail  phasing 
error  whan  the  actual  stub  length  being  slnulated  differs.  Effects 
howaver,  are  not  easily  dlscemable. 

6.  The  length  of  a atub  used  to  drive  the  bus  fron  a stub  location 
la  also  always  equal  to  a naln  bus  cable  elenent  In  length. 

7.  When  nultlple  stubs  are  slnulated,  all  stub  locations  being  used 
nust  contain  the  sane  nunber  of  nultlple  stubs.  For  exanple.  If 

10  stub  locations  arc  used  and  the  nultlple  stub  option  Is  elected, 
all  10  locations  would  have  the  sane  nunber  of  "M"  stubs.  If 
N ■ 4,  a total  of  40  stubs  would  load  the  bus. 

D.  COUPLER/STUB  CALIBRATION  DISCUSSION 


Coupler/Stub  characteristics  are  dctemlned  by  the  paranetcr  listed  under 
TRANSl  and  TRANS2  plus  ZL,  ZX,  REFl  and  assigned  line  length.  These 
parancters  describe  both  dynanic  and  steady  state  loading  effects. 

Dynanlc  loading  paraneters  consist  of  WNSQl,  TUOZVl,  ROOT2  and  the  capacitive 
load  Introduced  by  the  stub  length  LSTUB  array.  These  parancters  arc  dctemlned 
by  ncthods  dlscusitcd  In  Paragraph  II-D  and/or  with  aids  and  references  well 
known  to  control  systen  design  (eg..  References  10  and  11). 

Steady  State  loading  Is  dctemlned  by  REFl  (l.e.,  Fj^),  ZL,  GAINl  and  GAIN2. 

REFl  Is  calculated  fron  the  steady  state  Input  Impedance  presented  at  a 
stub  location  In  accordance  with  equation  14. 

GAINl  establishes  Che  voltage  levels  that  appear  on  the  stub  end.  This  para- 
ncter  setting  Is  not  generally  critical  and  can  be  made  equal  to  the 
steady  state  output  voltage  gain  of  the  transfomer  multiplied  by  0. 

Reflective  calculations  at  the  end  of  the  stub  adjust  for  resistive  loading. 

In  the  program,  the  parameter  REF4  is  calculated  from  ZL  using  equation 
14.  No  specific  entry  other  than  ZL  is  required. 


no  load  _ 
/|l+Rcf  l| 


GAIN2  establishes  th'*  steady  state  voltage  fed  back  onto  the  bus.  This  Is 
an  Important  Input  for  accurate  slmlation  with  numerous  stubs.  GAIN2  esn 
be  related  to  the  Insertion  loss  of  the  coupling  network.  An  equation  that 
can  be  used  to  provide  this  imput  Is  as  follows 


GAIN2 


oC  - (1  ♦ REFl) 
(1+REFl)  (GAINl) (F4) 


where  aC 
REFl 
GAINl 

F4 


- Insertion  loss  of  the  coupling  (voltage  ratio) 

- Input  parameter 

- Input  parameter 

- reflection  coefficient  at  end  of  stub  deteralned  via 
equation  14. 


aC  Is  a measured  Insertion  loss  of  the  coupler  and  equal  to  the  voltage  on  the 
bus  with  a loaded  coupler  divided  by  the  voltage  on  the  bus  with  no  coupler. 
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Multiple  stubs  are  simulated  by  providing  entries  that  duplicate  their 
loading  effects.  For  example,  REFl  is  entered  as  the  value  of  Fi  for 
the  number  of  parallel  couplers.  The  GAINl  entry  is  the  same  as  for  one 
stub.  The  GAIN2  entry  can  be  determined  from  the  equation  above  with 
adjusted  values  of  aC  and  REFl.  Parameters  effecting  dynamic  loading  are 
the  same  as  for  a single  stub. 
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V.  LABORATORY  EVALUATION 


This  section  describes  th«  bresdbosrd  hsrdvsre  used  for  sloulstor  verlfl- 
cstlon  end  provides  test  dsts. 

A.  DATA  BUS  BREADBOARD 

The  verlflcstlon  breadboard  Is  designed  In  accordance  with  MIL-STD-1553 
design  criteria.  It  consists  of  coupler  units,  cable  sections  of  various 
length,  connector  SMimted  resistor  terminators  and  a Bl-phase  L breadboard 
driver.  All  elements  contain  compatible  connectors  so  that  a variety  of 
data  bus  configurations  are  possible  without  soldering. 

Two  coupler  transformer  types  (l.e.,  two  different  transformers)  were 
purchased  by  the  statement  of  work.  Figure  27  provides  a coupler  schematic. 
Transformers  and  series  fault  Isolation  resistors  are  mounted  In  small  boxes, 
Pomona  model  2397  or  3230. 

The  transformers  were  arbitrarily  selected  to  have  markedly  different 
parameters.  The  transformers  and  key  characteristics  are  as  follows: 


PE  2244  (Pulse  Engineering) 

Turns  Ratio  1:2 

Primary  Inductance  6 MHy 

Prl/Sec  Leakage  Inductance  11  UHy* 

Winding  to  Winding  Capacitance  IS  pfd 

Primary  dc  Resistance  ^0  ohms 

Secondary  dc  Resistance  24  ohms 


PE  2226X  (Pulse  Engineering) 

Turns  Ratio  1:1 

Primary  Inductance  1.3  MHy 

Prl/Sec  Leakage  Inductance  1.5  UHy* 


* Measured  values . 
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winding  to  Winding  Cnpacltanco  3.6  pfd 
Prlanry  dc  Roolstanco  3.6  ohu 
Socondnrd  dc  Rcalntnnco  3.8  ohu 


Tranoforur  characterlatlca  art  aa  publlahcd  by  the  unufacturer  except  aa 
noted.  The  1:2  tranaforur  vaa  connected  to  provide  a atep  down  froa  the 
bua  to  a stub. 

The  cable  aectlona  are  fabricated  froa  68  oha  teflon  Inaulated  twisted 
shielded  pair  with  12  t%rlsts  per  foot.  Cable  Is  coapatlble  with  MIL-STD-1S53 
requlreunts  with  attenuation  characteristics  given  In  Figure  6.  Propo- 
getlon  delay  was  uasured  at  approxlutely  1.4  ns/ft.  Wire  la  ITT  type  2A 
SWTE  1936  AJTM.  Cable  sections  are  aa  follows:  eight,  10  feet  long;  seven, 

20  feet  long;  one,  30  feet  long;  one,  40  feet  long;  one,  120  feet  long;  and 
two,  2 Inches  long. 

Teralnatlon  resistors  Include  eight  2.2  Koha  resistors  for  reute  teralnal 
load  slaulatlon  and  two  68  oha  resistors  for  cable  end  teralnatlon. 

The  Bl-phase  L breadboard  driver  circuit  produces  three  level  signalling 
coapatlble  with  M1L-STD-1S53.  The  circuit  la  housed  In  a 6"x4"x2"  chassis 
with  an  output  chassis  connector  that  utes  with  the  connectors  provided  on 
uch  end  of  the  cable  sections.  Figure  28  provides  a scheutlc  Identifying 
coaponent  parts.  Figure  29  shows  Input  wavefora  required  for  operation.  The 
breadboard  driver  also  requires  external  power  supply  Inputs  of  <fS  volts, 

100  allllaaps  and  +8  to  'MS  volts  at  225  alllleaps. 

A laboratory  pattern  generator  Is  used  to  operate  the  driver.  DATAPULSE 
Models  203  and  208  arc  two  exaaples  of  coapatlble  pattern  generators. 

Inputs  to  the  breadboard  driver  arc  TTL  coapatlble.  Figure  29  defines  wevc- 
fora  Inputs  and  outputs  of  the  driver. 
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Figure  2*1.  Breadboard  Driver  Input  and  Output  Waveforms 


B.  TEST  DATA 


Basically  two  test  configurations  ifere  used  for  coaparing  breadboard  and 
siaulator  outputs.  The  first.  Setup  A (Figure  30),  was  used  for  characterizing 
transfozaers  and  stub  loading  effects  of  different  lengths.  The  second. 

Setup  B (Figure  31),  was  used  for  evaluating  eight  stub  configurations. 

In  Setup  B,  the  generator  was  swved  to  different  locations  which  included 
driving  the  bus  froa  stub  ends  as  is  nonsally  the  case  in  a NIL-STD-1553 
application. 

Figure  32  shorn  a unique  case  of  driving  120  feet  of  cable  terainated  in  68 
ohaw.  The  figure  shows  photographs  taken  of  wavefonss  as  they  appear  at  both 
cable  ends  along  with  siaulator  outputs  froa  the  identical  configuration. 

Minor  distortion  noted  at  the  breadboard  generator  is  believed  caused  by 
cable  characteristic  lapedance  variation  with  frequency.  End  termination 
coaparlsons  agree  favorably.  The  exponential  rise  effect  illustrates  skin 
effect  and  its  slaulation. 

Figures  33  to  38  reflect  data  from  the  coupler  that  uses  2:1  transformers 
(PE  2244).  Figures  33  *nd  34  show  Setup  A with  different  stub  lengths  and 
consequently  the  effects  of  the  capacitance  stub  loading.  Figures  35  and 
36  provide  comparison  ffor  eight  stub  loads  of  Setup  B with  generator 
driving  from  the  left  end.  Figures  37  and  38  are  similar  cases  except 

driving  from  Stub  #1,  forty  feet  from  the  left  end  of  the  bus. 

Figures  39  and  40  show  data  compar'*8ons  for  the  1:1  transformer  (PE2228X) 
coupler  units.  Figure  39  is  a one  20  feet  case  of  Setup  A.  Figure  40  is 
the  eight  stub  configuration  of  Setup  B. 

Figure  40  indicates  that  paraaeter  or  even  a model  change  would  improve  simu- 
lation of  data  froa  generator  driving  from  stub  to  bus. 
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Figure  31.  Setup  B 
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rifure  32.  Sl«iUtor/Lab  C<»parlaon  for  120  feet  oi  cable  Into  6S 


Figure  33.  Slnulator/Lab  Conparlson  for  one  10  foot  stub 
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Flgura  34.  Slaulator/Lab  Coaparlson  tor  ona  20  faot  atub 
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Figure  35  • Slaulator/Leb  Coapariaon  for  eight  10  foot  atuba 


Slailacor 
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Figure  36.  Slaulator/Leb  Coaparlaon  for  eight  10  foot  acuba 


Figure  37.  SlmuXator/Lab  Comparison  driving  Stub  #1 


SlMlator 


Figur*  38.  St*ul«tor/L«b  Coaparison  at  Stuba  #1  and  #8 


Figure  39.  Slaulator  Lab  Coaparlson  for  one  20  foot  atub 


Breadboar«< 
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Figure  40.  SlaulaCor/Lab  Coaparleon  for  eight  10  foot  atube 


VI. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  network  simulator  developed  on  this  program  should  prove  a valuable  tool 
for  MIL-STD-1553  application  study  and  for  designing  couplers.  Furthermore, 
since  the  program  Is  FORTRAN  coded,  modifications  to  the  basic  software  could 
and  should  broaden  Its  use. 

Because  the  contract  scope  was  limited,  several  areas  were  not  investigated 
as  fully  as  desired.  These  Include  the  characterization  of  transformers  and 
stubs.  An  original  goal  was  to  fully  simulate  transformers  from  resistance, 
inductance  and  capacitance  measurements  or  via  similar  manufacturer's  specifica- 
tion data.  As  Indicated  In  Paragraph  II-D,  this  goal  was  not  achieved.  Addi- 
tional work  in  this  area  Is  recommended.  Additional  stub  analysis  is  also 
recomnended . 

The  transmission  line  algorithm  that  takes  skin  effect  into  account  is  believed 
a noteworth  development.  With  three  parameters  obtainable  from  laboratory 
data,  an  accurate  line  simulation  Is  achievable  for  a wide  range  of  transmission 
lines.  Simulation  accuracy  can  be  further  improved  by  adding  terms  and/or  by 
changing  calibration  points. 
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